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Abstract. Results of identified charged hadrons produced in Au+Au collisions at 
yfsm = 130 GeV as a function of centrality are discussed. The (pj) of protons and 
antiprotons is significantly larger than in p-p collisions at the same center of mass en- 
' er gy- Its scaling with the number of participant nucleons, an increase of 20 ± 5% from 

peripheral to central, may be due to radial expansion. In a hydrodynamical picture a 
simultaneous fit of pions, kaons and protons is consistent with a radial flow velocity 
^ i , of Pisurface = 0.72 ±0.01 and a temperature of T = 122 ±4 MeV for the 5% most 

{S) ' central collisions. Dominant proton yields at pj > 1 .5 GeV/c and antiproton yields at 

, pj > 1.8 GeV/c may be explained by a combination of suppression of pion yields at 

' high pj and the boost of low pj protons to high pj from radial flow. 
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1. Introduction 



It is predicted by theoretical calculations that highly excited and dense hadronic matter 
undergoes a phase transition from a hadronic gas to a deconfined state of quarks and gluons 
called a Quark Gluon Plasma (QGP) [ [l]]. A few ps after the Big Bang [ 0], before nuclei 
and atoms were formed, the very hot and baryon free early universe may have existed in 
such a state of deconfinement. 

The objective of relativistic heavy-ion physics is to produce an energy density in the 
overlap region between two colliding ions that is high enough to probe the phase transition 
between hadronic matter and quark matter. The quarks and gluons coexisting in a QGP state 
cannot be measured directly. The experimental challenge is to study a suite of observables 
that are indirect probes of plasma formation. The one observable that is measured and 
discussed here are the transverse momentum (pj) spectra of produced hadrons, namely: 
pions, kaons and (anti)protons. Unlike photons and dileptons, the hadrons decouple late in 
the reaction due to strong interactions; therefore, the momentum spectra are sensitive to the 
system at freezeout. 
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If the collision can be described as an expanding, relativistic fluid, then hydrodynamic 
models with various equations of state can be used to predict the hadron spectra [ph. If both 
hydrodynamic behavior and longitudinal boost invariance are assumed, then the amount of 
expansion generated in the collision can be extracted from the data at low momenta (pr < 
1 .0 Ge V/c) using a hydrodynamics-inspired parameterization. Previously, the expansion 
of matter in Pb+Pb collisions at CERN SPS energies was similarly obtained from hadron 
spectra and HBT data [ Q Bfl. In addition, the measured yield of high pj hadrons may 
reveal the amount of partonic energy loss when compared to calculations in perturbative 
Quantum Chromodynamics (pQCD) parton models [ ^, M], 

The hadrons from Au+Au collisions at a roughly ten times higher energy density are 
measured at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory 
(BNL). The PHENIX experiment at RHIC is described briefly in Section f| 

In the first year of data taking, the identified charged hadrons produced in five differ- 
ent event centralities at the collision energy of 130 GeV in Au+Au collisions are measured: 
0-5%, 5-15%, 15-30%, 30-60% and 60-92% of the total inelastic cross section of 6.8 b. 
The data reduction, corrections, extrapolations, and systematic uncertainties are summa- 
rized in Section |[ The centrality selected spectra, yields, mean transverse momenta ((pr)), 
and total charged multiplicity are presented in Section |[ An interpretation of the data based 
on radial, hydrodynamic expansion is provided in Section |[ The implications of the results 
are discussed in Section ^. 

2. Experiment 

2.1. The PHENIX detector 

The PHENIX [ [| ^ detector is designed for the measurement of e, p, y, and identified 
hadrons over a large momentum range for pj > 0.2 GeV/c. It is comprised of four spec- 
trometer arms: two at mid-rapidity called the central arms ("east" and "west") and two at 
more forward rapidities called the muon arms ("north" and "south"). The central arm ac- 
ceptance is |t]| < 0.35 with each arm subtending Ji/2 in azimuthal angle ((j)). In the first year 
of physics running in 2000 (Run-1), half of the central arm detectors were instrumented, as 
shown in Fig. |[ 

The two types of collision axis detectors, the Zero Degree Calorimeters (ZDC) and 
the Beam-Beam Counters (BBC) measure both the initial time and vertex position of each 
event. The ZDC resolution was Oe/E < 20% at the single neutron energy of E n = 100 GeV 
with an acceptance of |t|| > 6. The BBC resolution was a z = 1.5 cm along the beam axis 
and o, = 70 ps in timing, with an acceptance of 3.0 < r| < 3.9. 

The time-of-flight detector (TOF) in the east arm subtends 7t/4 in <j) and was used 
for the identification of 7C , K*, and (anti)protons. The TOF measures the arrival time of 
charged hadrons at the TOF wall positioned 510 cm from the collision vertex. In Run-1, a 
timing resolution of 1 15 ps was obtained. 

In Run-1, PHENIX had two layers of multi-wire proportional chambers with pad- 
readout (PCI, PC3) that measure a three-dimensional space point of a charged track. The 
points from PCI are used in the global track reconstruction to determine the polar angle 
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PHENIX Detector - First Year Physics Run 




West BeamView 



Fig. 1. A cross-sectional view of the PHENIX detector, transverse to the beam-line. Within 
the two central arm spectrometers the detectors that were instrumented and operational 
during the y/sm = 130 GeV run are shown. 

of the track. The PC3 positioned in the east arm is used to minimize the background 
contribution of albedo and non-vertex decay particles; it is used in the high-p^ charged 
hadron analysis [ |7|. 

The drift chambers measure the track curvature between 202 and 246 cm radially, for 
which the particle momentum is determined. In Run-1, a resolution of o p / p = 0.6% ©3.6% 
was obtained. The absolute momentum scale was known better than 2%. 

3. Data Reduction and Analysis 

3.1. Event sample 

A total of 5M minimum bias triggers were recorded at ^/snn = 130 GeV in the PHENIX 
ZDC detectors. The event centrality was determined using a correlation measurement 
between the neutral energy deposited in the ZDCs and fast charged particles (p > 0.69) 
recorded in the BBCs [ O], In order to minimize the amount of albedo from the central 
magnet pole faces, positioned at ±40 cm, only those events with a vertex distribution within 
±30 cm were selected. Additional event criteria include the magnetic field setting and de- 
tector stability. The resulting event sample, 140000 minimum bias events, was 92 ±4% 
of the total inelastic cross section of 6.8 b [ [TT| ]. The events were subdivided into five 
centrality classes: 0-5%, 5-15%, 15-30%, 30-60% and 60-92% of the total geometric 
cross section. In each centrality selection, the number of participating nucleons (N part ) and 
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number of binary nucleon-nucleon collisions (A^ co n) were determined using a Monte Carlo 
Glauber model calculation [ O]. 

3.2. Particle identification 

Particle identification for charged hadrons was performed by combining the information 
from the tracking detectors with the timing information from the BBC and the TOE Tracks 
at 1 GeV/c in momentum point to the TOF with a projected resolution o pro j of 5 mrad in 
azimuthal angle and 2 cm along the beam axis. Tracks that point to the TOF with less than 
2.0 O pro j were selected. 

The measured momentum (p), path length (L), and time of flight (t ) in the spectrometer 
were used to calculate the particle mass, with p = L/ct in Eq. ([!]). 




The width of the mass-squared distribution o m i was calculated as a function of momentum, 
using the known detector resolutions and the relation in Eq. (Q), 

o\ = C, 2 • 4m 4 (l + ^ + C 2 2 • 4^ + C 3 2 • [4/(m 2 + p 2 )} . (2) 

The constants are related to the measured drift chamber momentum resolution (C\, Cj) and 
combined flight path length and flight time errors (C3). Tracks that fell within 2a m i of a 
measured mass centroid were identified. 

3.3. Corrections 

The raw spectra include inefficiencies from detector acceptance, resolution, particle de- 
cays in flight and track reconstruction. Corrections for these effects were determined using 
a Monte Carlo simulation of the detector response. The baseline efficiencies were de- 
termined by simulating and reconstructing single hadrons. Multiplicity dependent effects 
were evaluated by embedding simulated single hadrons into real events, and measuring the 
degradation of the reconstruction efficiency. For each hadron, sufficient statistics were gen- 
erated fiat in azimuthal angle and rapidity for \y\ < 0.5 for every measured pr-bin such that 
the statistical error in the data dominates. The statistical errors from the correction factors 
were added in quadrature to the statistical errors in the data. Corrections for feed-down 
from weak decays were not applied, but an upper limit is estimated based on a MC simu- 
lation of A particles within the PHENIX acceptance. More details on the corrections are 
reported elsewhere [ [l3j [l4j]. 

3.4. Extrapolations 

The yield integrated over transverse momentum, dN/dy, and the integrated average trans- 
verse momentum, (pr), were determined for each particle after fitting and integrating func- 
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tions that describe the spectral shape. A power-law and mj exponential were fit to the pion 
data, while pj and 017- exponentials were fit to each kaon and proton spectrum. 

3.5. Systematic uncertainties 

The overall uncertainty on dN/dy is 13%, 15% and 14% for pions, kaons and protons, 
respectively. Uncertainties on (p T ) depend on the extrapolation and background uncertain- 
ties; the uncertainties are 7%, 10% and 8% for pions, kaons and protons, respectively. An 
estimated 24% of the measured proton and anti-proton yields are due to A weak decays, 
based on a Monte Carlo simulation. 



4. Results 



4.1. Transverse momentum distributions 



The 7t , K , p and p invariant yields as a function of transverse momentum (pr) for the 
most central, mid-central, and the most peripheral events are shown in Fig. ^. These data 
were reported previously [ OJ . 



10 r 
10 r 
1 r 
10 r 

-2 : 
10 r 

-3 : 

10 i 

in 

10 I 
1 

10" 1 i- 

10 
-in 

10" k 

10 

1 

10 1 
10 2 
10- 3 t 

10 



.""t^ positive 0-5% j. n cj, negative 0-5% 



• P* 



^8 



15-30% 




60-92 % 



□ E. 4 I. 



15-30% 



* HP 



"4 □ 

1 



8 □ 



Fiig. 2. The invariant yields for pions, 
60-92 % k aonSi an( j (anti)protons produced in 
three different classes of event cen- 
trality: 0-5% (top), 15-30% (middle) 



0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 

p T (GeV/c) p T (GeV/c) 



I * t an d 6 °- 92% (bottom) [[l3|]. 



Transverse momentum spectra of non-identified charged hadrons and neutral pions 
were reported by PHENIX in the range 1 < p T <5 GeV/c [g]. In Fig. |, the 71° data are 
shown with the charged pion data for the most central events. The comparison confirms the 
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apparent baryon to meson dominance at pj w 1.5 GeV/c and pj sw 1.8 GeV/c for positive 
and negative particles, respectively. 
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Fig. 3. The invariant yields as a function of pj for the positive (left) and negative (right) pi- 
ons in the most central events are consistent with the Jt° measurement from an independent 
PHENIX analysis, verifying the apparent dominance of baryon yields at high pj. 



4.2. Yield and (pj) 

The dN/dy is normalized to pairs of participant nucleons for each centrality selection, as 
shown in Fig. [|. The dashed lines are the systematic uncertainties in the Glauber model 
calculation of the number of participant nucleons. The kaons and (anti)protons are scaled 
for clarity. The error bars include both statistical and systematic errors. 

The mean transverse momentum increases with the number of participant nucleons 
by 20 ±5% for pions and protons, as shown in Fig. 0. The (p?) of particles produced 
in p-p and p-p collisions, interpolated to RHIC energies, are consistent with the most 
peripheral pion and kaon data; however, the (pj) of protons produced in Au+Au collisions 
is significantly higher. This dependence on the number of participant nucleons may be due 
to radial expansion. 

4.3. Charged particle multiplicity 

The measured dN/dy for each hadron species is converted to dN/dr\, and the total dN/dr\ 
is calculated after summation. In Fig. 0, the resulting dN /dry per participant nucleonpair is 
within 5% when compared to the measurement made by PHENIX using the pad chambers 
alone [ 111]]. The results from other RHIC experiments are also shown as a comparison. 
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Fig. 4. The integrated dN/dy per participant pair for K + , K + , p (left) and jr - , K - , and p [ 
|j~3| ]. The dashed lines (left) are the systematic uncertainties in the Glauber model calcula- 
tion Of JVpart. 




50 100150 200 250 300 350 

Npar, 



50 100150 200 250 300 350 

Npar, 



Fig. 5. The integrated mean pj for pions, kaons, and (anti)protons produced in the five 
different classes of event centrality [|l3j]. The error bars are statistical only. The systematic 
uncertainties are 7%, 10% and 8% for pions, kaons and (anti)protons, respectively. The 
open points are equivalent average transverse momenta from p-p and p-p data, interpolated 
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Fig. 6. The total charged multiplicity (open) in Refs [ |17|, [18], 19 1 and the total identified 
charged multiplicity (closed) scaled by the number of participant pairs are plotted as a 
function of the number of participants. 



5.1. Parameterization 

Radial flow would impart a velocity boost to the hadrons. Heavy particles should be 
boosted to higher pj, depleting the cross section at lower pj, and resulting in a higher 
inverse slope in the spectrum. In order to measure the amount of radial expansion, a hydro- 
dynamics inspired parameterization is used to fit all spectra simultaneously. The functional 
form used describes a boosted thermal source, and is based on relativistic hydrodynamics, 
see U. Heinz et al. in [[T^]. The particles are assumed to collectively expand with a velocity 
profile that is linear with the radial position. The particle density distribution is assumed to 
be a Gaussian in the radial position. In addition to longitudinal boost invariance, it is also 
assumed that all particles decouple kinematically on the freezeout hypersurface at the same 
freezeout temperature Tf a . The transverse velocity profile is parameterized as: 

Pr(§) = P^", (3) 

where t, — r/R, p s is the surface velocity, n is the profile, and R is the radius of the ex- 
panding source at freezeout (0 < Z, < 1 in the case of a flat particle density distribution) [ 
jig ]. Each fluid element is locally thermalized and is given a transverse rapidity boost p 
that depends on the radial position as: 



p=^tanh- 1 ((3 i ) 



(4) 
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The mj dependence of the yield dN jmjdmj is calculated after integrating over the az- 
imuthal, rapidity, and a final numerical integration over the radial components of the source 

£L- r A J mT f { m (^^)/o(^^)^, (5) 

using the modified Bessel functions Iq and K\ . The parameters in Eq. (|J) are the freezeout 
temperature 7f , the normalization A and the maximum surface velocity p s . A linear veloc- 
ity profile (n = 1) within the source is used. This choice was guided by the profile observed 
in a full hydrodynamical calculation by Kolb and Heinz.* 

5.2. Simultaneous fits to data 

The amount of radial expansion is determined by fitting all particle transverse momentum 
spectra simultaneously using Eq. for pj < 1 .0 GeV/c. We make a grid of combinations 
of temperature and velocity, and perform a chi-squared minimization to extract the normal- 
ization, A, for each particle type. The fit is done simultaneously for all particles, with % 2 
contours that show the anti-correlation of temperature and velocity. The resonance region 
in the pion spectra is excluded in the fit,* as was done in S+S collisions at the CERN SPS 
[ [Z0]]. The resulting fits are plotted with the spectra for all centralities in Fig. [|. 

We find that the transverse momentum distributions are consistent with an expanding 
system with a surface velocity of p s — 0.72 ± 0.01 which decouple at a common tempera- 
ture of T fo = 122 ±4 Me V. For the 5-15% centrality, the fit results are r fo = 125 ± 2 MeV 
and p\ = 0.69 ±0.01. 

In [ 0], the hadron multiplicities at high pr are different for mesons and baryons. The 
suppression of pion yields at high pj and the boost of low pj protons to high pj from 
radial flow may explain the dominant proton and antiproton yields at pj > 1.5 GeV/c and 
1.8 GeV/c, respectively. 

5.3. Systematic uncertainties 

We find that Eq. (|) is consistent with a full hydrodynamic calculation when a linear ve- 
locity profile is used [ Q. Nevertheless, a parabolic profile results in an increase of «13% 
in p v and «5% in Tf Q . A constant particle density distribution used with a linear velocity 
profile decreases p s by «2% with a negligible difference in the temperature Tf Q . As a test 
of the assumption that all the particles freezeout at a common temperature, the simulta- 
neous fits were repeated for pions and protons only. The difference in 7f is negligible 
within the measured uncertainties. The effect of including the pion resonance region in the 
simultaneous fit decreases Tf by 20 MeV. 

6. Conclusions 

Identified charged hadrons produced in Au+Au collisions at 130 GeV in five different cen- 
trality selections have been measured. Proton and antiproton yields at high pj dominate 
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Fig. 7. The hydro-inspired parameterization is fit simultaneously to the data in each event 
centrality 



over pion yields in central collisions when compared to the PHENIX measured Jt data up 
to pt — 4 GeV/c. A combination of suppression of pion yields at high pj and the boost of 
low p T protons to high p T from radial flow may be an explanation. 
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